Non-local composite spin-lattice polarons in high temperature superconductors 
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The non-local nature of the polaron formation in t-t'-t"-J model is studied in large lattices up to 
64 sites by developing a new numerical method. We show that the effect of longer-range hoppings 
t' and t" is a large anisotropy of the electron-phonon interaction (EPI) leading to a completely 
different influence of EPI on the nodal and antinodal points in agreement with the experiments. 
Furthermore, nonlocal EPI preserves polaron's quantum motion, which destroys the antiferromag- 
netic order effectively, even at strong coupling regime, although the quasi-particle weight in angle- 
resolved-photoemission spectroscopy is strongly suppressed. 

PACS numbers: 71.10.Fd, 71.38.-k, 02.70.Ss, 75.50.Ee 



It is a matter of long debates whether electron-phonon 
interaction (EPI) plays an essential role in the forma- 
tion of the puzzling properties of high temperature su- 
perconductors Early theoretical studies of the an- 
gle resolved photoemission spectroscopy (ARPES) of un- 
doped cuprates were based on the t-t'-t"-J model, where 
doped holes into an antiferromagnetic (AFM) material 
are characterized by exchange constant J and hoppings 
up to 3rd near neighbors(NN) with amplitudes t, t' and 
t". This approach successfully described the dispersion 
of the experimentally observed peak Q and appeared to 
exclude strong EPI. However, there was a strong contra- 
diction between the theoretical and experimental line- 
shapes: in contrast with the very broad peak in exper- 
iments flj the theory predicted a sharp peak [1, 0, @1- 
This contradiction was resolved under assumption that 
the undoped parent compounds are in the strong cou- 
pling regime (SCR) of the EPI. In the simplest case of t- 
J-Holstein model, where a hole interacts with dispersion- 
less optical phonons through on-site local coupling, it was 
shown that at SCR the spectral weight of the quasiparti- 
clc is almost completely transferred to the broad Franck- 
Condon shake-off peak which inherits the dispersion of 
the hole when it does not interact with phonons Fur- 
thermore, the inheritance of the dispersion of the nonin- 
teracting quasiparticle by the Franck-Condon peak was 
shown to be a property of a broad class of models with 
EPI Q • Recently, the relevance of the t- J-Holstein model 
for cuprates has been rather convincingly questioned Q. 
It has been shown that SCR leads to a picture of a local- 
ized static hole with 4 broken bonds about it. In this case 
the percolative model predicts that AFM order survives 
up to critical concentration x w 0.5 in contradiction with 
the experimental value x ~ 0.02 — 0.04. Another problem 
of the t-J-Holstein model is the very large effective mass 
in the whole SCR 0] contradicting transport and optical 
properties of lightly doped cuprates 0] . 

The above arguments, in favor or against the SCR in 



undoped cuprates, suggest that the answer must be found 
in a class of more realistic t-t'-t"-J models where EPI, as 
it is in cuprates is nonlocal. Besides of the particu- 
lar relevance of the above models in the physics of high 
temperature superconductivity, there is a general inquiry 
asking whether long-range interaction can introduce pro- 
found qualitative difference compared with t-J-Holstein 
model. However, study of the t-t'-t"-J model with nonlo- 
cal EPI by exact diagonalization (ED) and Diagrammatic 
Monte Carlo (DMC) methods is hindered by their lim- 
itations. ED method has been restricted to only small 
two-dimensional systems 0, [lH and the nonlocal effects 
have not been studied so far. On the other hand the study 
of the incoherent motion of the hole requires transforma- 
tion of the DMC method [6| into the direct space where 
DMC encounters the sign problem Q. 

In this Letter we app ly recently introduced coherent 
state basis approach [l2| to the t-t'-t"-J model with non- 
local EPI, and develop a Coherent States Lanczos (CSL) 
method which is able to provide a reliable description 
of the ground state properties of two-dimensional sys- 
tems up to 64 sites. We validate it by successfully com- 
paring the results with exact DMC data. We find that 
in the t-t'-t"-J model with local Holstein EPI, in con- 
trast with t-J-Holstein model, the influence of the EPI 
on the properties of nodal k = (7r/2,7r/2) and antinodal 
k = (tt, 0) points is significantly different introducing a 
considerable anisotropy driven by long-range hole hop- 
pings: the antinodal states are in the SCR with nearly 
zero quasiparticle weight whereas the quasiparticles in 
the nodal point are only weakly dressed by phonons in 
a very broad range of EPI. We also show that nonlocal 
EPI qualitatively changes the basic features of the com- 
posite spin-lattice polarons. First, according to general 
result [l3|, effective mass in the SCR of nonlocal EPI 
is considerably lighter than that of a model with local 
EPI. Second, the motion of the hole over the NN is not 
strongly suppressed by nonlocal EPI and, thus, the ten- 
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dency toward robust AFM order up to high doping levels, 
found in t-J-Holstein model, does not survive in the more 
realistic t-t'-t"-J model with nonlocal EPI. 

The minimal Hamiltonian for t-t'-t"-J model with ex- 
tended EPI is a sum of t-t'-t"-J Hamiltonian H^^/^'i j 
@i S, [3 1 Hamiltonian Hh-ph, and Hamiltonian of 
dispersionless phonons Hph = uiq J2i c\ci (cj is the cre- 
ation operator of a phonon at site i with frequency loq). 
In the spin wave approximation the t-t'-t"- J Hamiltonian 
reads 



tt't" J 



Y.Y.[Y. Ak,jlh,+sa,- + Ah,hlf,+sbir + h.c. j , 
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where A (B) is the sublatticc with spin up (down) and 
at (6t) is the operator creating magnon in the A (B) 
sublatticc with dispersion Qq 0, 0, Q • The opera- 
tor in the direct space fj {hi) creates a spinless hole 
on the site i of the sublatticc A (B) . The bare hole 
dispersion is = 4t cos((7a;) cos(gy) -|- 2t {cos{2qx) + 
cos(2(7y)) and the hole-magnon coupling is Mq^i ~ 
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ti^e"?'^'+^ ) , where itq- and are the 



Bogoliubov factors, S is an unitary vector connecting NN, 
and N is the number of lattice sites. The sum over q is 
restricted inside the magnetic Brillouin zone. The EPI 
Hamiltonian 
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in defined in terms of Holstcin. g(0) ~ g, and nonlocal 
coupling to the NN lattice displacements g{5) = gi. Be- 
low we set h = 1 and t ~ 1. 

In order to study the ground state of this model, we use 
CSL procedure based on the Lanczos recursion method 
that, starting from the hole in the quantum Neel state 
without excited magnons, lO)*"', generates the subspace 
spanned by the basis 



|j,/il,...,/iAr,gi,...,(f/,0 = \h) 
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where \h)j indicates the state with the hole on the 
site j, and i runs over the lattice sites. The nota- 
tion \qi, ....,qi,l) labels the standard Bonfim-Reiter l- 
magnon states which arc sufhcient to reproduce the re- 
sults of self-consistent Born approximation (SCB A) [l^ . 
These states are given in terms of an ordered product 
of at and 6t operators acting on the magnon vacuum 

|0)<"' (kl,....,<7l,0 = 
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FIG. 1: (Color online) Comparison of energy (a), spectral 
weight (b), phonon (c) and magnon (d) average numbers 
(squares) obtained by CSL method at A'' = 64 with those 
obtained by DMC method (line) for the t-J-Holstein model 
at J/t = 0.3, Luo/t = 0.1, and k = (7r/2,7r/2). 



&tat ....lO)'"' ) if the hole is on A (B) sublatticc. 
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Within the Bonfim-Reiter basis, the Lanczos or any other 
ED procedure is capable of reproducing the ground state 
properties and spectral functions within the SCBA if 
the following approximations are implemented: 1) the 
magnons involved into the string at6t_^....|0)<"' are 
considered to be distinguishable, leading thus to the 
unity normalization of the state gj, 2) the hole- 

magnon coupling term, as in the retraceable-path ap- 
proximation [lij . is considered to be able to annihilate 
only the last excited magnon. For the study of the 
ground state properties the situation is considerably sim- 
plified because for J/t > 0.3 the states containing less 
than 4 magnons are sufficient to reproduce the results 
of SCBA [15[. The real bottleneck for increasing the 
size of the system comes from the phonon basis states 



|0),^''' , where |0),-''' is the phonon vacuum 



at the site i. The number of relevant phonon states 
quickly grows at intermediate values of the EPI strength 
restricting the lattice sizes to 10 sites [lH, the maximal 
size we are aware of, for quantum phonons and to 20 sites 
in the fully adiabatic limit Q. To circumvent the bot- 
tleneck we switch from the basis of many-phonon states 
to that expressed in terms of coherent states (CS) 
[12| which are nothing but the canonical transformations 
for independent oscillator model 



\h,t) 



^gh(b,-b\) jgyiph) 



?i=0 



(2) 



They are renormalized by different parameters h. For 
h = Q the CS is the bare state without EPI and for 
/i = 1 it is the exact solution of the independent oscilla- 
tor model with local coupling g. In principle, this substi- 
tution does not introduce any truncation in the Hilbcrt 
space since, varying h in the complex plane, the local ba- 
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FIG. 2: The spectral weight (a and c) and the mean kinetic 
energy (b and d) as function of A for nodal (circles) and antin- 
odal (triangles) points. The parameter values are: J /t = 0.4, 
uo/t = 0.2, gi=Q, N = 64. 



FIG. 3: (Color online) Average number of phonons (a and 
c) for nodal (circles) and antinodal (triangles) points; PPDF 
(b and d) for nodal (solid lines) and antinodal (dashed lines) 
points. The parameter values are: J/t — 0.4, uo/t = 0.2, 
gi=0, N = 64. 



sis ([5]) is ovcr-complctc. Naturally, particular realization 
of the CSL approach requires a truncation. Actually only 
few CS arc enough to reproduce the ground state prop- 
erties of the composite spin-lattice polaron. Comparison 
with exact DMC data shows that the following truncation 
scheme minimizes the computational efforts with high 
accuracy for results. A finite number Mph of CS is cho- 
sen, characterized by real values ha {a = 0, ...,Mph — 1) 
equidistant in the range [0, 1], and only the phonon states 

Y[f=i l^^ai, i) with X^iLi — ^'Ic,ph arc included into the 
basis. Besides, taking advantage of the fact that the 
hopping and EPI are restricted to 3rd NN and NN re- 
spectively, we restrict the states with nonzero parame- 
ters h to 3rd NN from the hole. Comparison (Fig. 1) 
of CSL results for t-J-Holstein model on 8 x 8 lattice 
with approximation-free result by DMC in the thermo- 
dynamic limit shows that our method at M^ h = 4 
and Mc^ph = 3 is valid for all ranges of EPI [17|. The 
small differences in the weak EPI regime are mainly due 
to finite size effects. At large dimcnsionless couplings 
A = g^LO^/At the polaron size becomes small and the 
agreement becomes excellent. Our method perfectly re- 
produces crossover between weak and strong coupling 
regimes at critical value Ac ~ 0.4 where i) the spectral 
weight goes to zero, indicating suppression of the coher- 
ent motion; ii) the lattice distortions rapidly increase; 
and iii) the size of magnetic polaron reduces. 

First, we extend t-J-Holstein model with local EPI 
adding next NN, t', and next next NN, t", hoppings. 
Comparing physical properties of these models we ob- 
serve that the longer-range hoppings introduce a strong 
anisotropy. Physical properties of the states in the nodal 
(7r/2,7r/2) and antinodal (tt, 0) points are very similar 
in the t-J-Holstein model but considerably different in 
the tt't"J-IIolstcin model (Fig. 2 and 3), as it is evident 
comparing Z- factor (Fig. 2a and 2c), mean kinetic energy 



(Fig. 2b and 2d), and average number of phonons (Fig. 3a 
and 3c) . There is a wide range of the EPI strengths from 
the weak to intermediate coupling regime where the hole 
is a heavily dressed polaron in the antinodal point but an 
almost free particle in the nodal one. Comparing spectral 
weights and kinetic energies of the two models we con- 
clude that the longer-range hoppings reinforce coherent 
motion in the nodal point and suppress it in the antinodal 
one. To characterize the anisotropy induced by t' and t" 
we calculate the phonon probability distributions func- 
tion (PPDF) P{x) = {^GS \x){x\ ^Gs), where *gs is 
the wave function of the ground state and x is the lat- 
tice distortion of the site where the hole is, just below the 
crossover into the SCR (Fig. 3b and 3d). We find that for 
the t-J-Holstein model both nodal and antinodal points 
display a bimodal structure characteristic of the inter- 
mediate coupling regime, where the ground state of the 
system is a quantum mixture of localized and delocal- 
ized states [18[. To the contrary, in the tt't"J-Holstein 
model the antinodal point has bimodal PPDF while that 
at the nodal point displays the maximum at small lat- 
tice distortions x, characteristic of weak coupling regime. 
The resulting picture of the anisotropic EPI, consider- 
ably stronger near the antinodal point than at the nodal 
one, is expected to be qualitatively similar in the case 
of an underdoped system and is in line with experimen- 
tal findings. The explanation of the vertical dispersion 
of the ARPES near the antinodal points requires large 
Z (weaker coupling) along the nodal direction and much 
smaller Z (stronger coupling) along the antinodal ones 

El. 

To study the role of nonlocal EPI we compare the 
properties of the tt't"J-Holstein model with that includ- 
ing nonlocal coupling to the NN lattice displacements 
gi =5/2 (Eq.([T|)). In both cases the decrease of the Z- 
factor (Fig. 4a) indicates the suppression of hole ground 




FIG. 4: (Color online) The spectral weight (a), the spin deviation around the hole (see text) (b), and the absolute value of the 
average kinetic energy for magnon assisted NN hoppings, (c), for gi = (diamonds) and gi — 0.5g (crosses). The parameter 
values are: J/t = 0.4, ujo/t = 0.2, t = —0.5t, t — OAt, N = 16. Critical Ac is defined as that where the spectral weight is 
reduced by 70% with respect to that at A = 0. 



state contribution to ARPES. However, the increase of 
the EPI range slightly enhances (reduces) the suppres- 
sion of the Z-factor at A < Ac (A > Ac). On the contrary, 
when EPI is non local the polaron effective mass, that 
is a measure of the coherent motion, becomes larger in 
the weak coupling regime and smaller in the SCR. In the 
latter case, the mass for nonlocal EPI at A = l.lAc is one 
order of magnitude lighter: diagonal nid (along = ky) 
and transverse rrit (along k^ = —ky) masses for local EPI 
are = 71.4 and mt = 65.7 whereas for nonlocal EPI 
their values are irid = 7.1 and mt = 7.5, respectively, in 
units of the bare effective masses (g = gi = 0). Simi- 
lar effects have been discussed in a number of different 
models 

Figure 4b and 4c show the dependence of spin devia- 
tion and the contribution to the kinetic energy arising 
due to magnon assisted NN hoppings, Kt, on the di- 
mensionless coupling constant A. Spin deviation, SD = 
(Safm — {Snn))/Safm, measures how much the neigh- 
boring spin to the hole, Snn-, deviates from its value in an 
ideal antiferromagnet (Safm)- Absolute value of i^^t and 
SD are the measure of the intensity of the NN hopping 
and include both coherent and incoherent contributions. 
The decrease of the two above quantities, (Fig. 4b and 
4c) , signals on the suppression of the motion over NN as 
expected in presence of holc-phonon interaction. How- 
ever, by increasing the EPI range, the motion of the hole 
on NN is suppressed by the nonlocal EPI much less ef- 
fectively than it is in the case of the Holstein coupling. 
This effect is one of the main results of this work. It 
can be easily interpreted in the SCR where the notion of 
the adiabatic potential makes sense. For the Holstein lo- 
cal coupling the self-consistent adiabatic potential for the 
hole is a deep J-function, preventing, thus, the hole from 
the motion over the NN. To the contrary, the motion of 
the hole in case of a long-range EPI is not restricted to 
a single site, preserving, thus, the motion over the NN 
even at SCR. This allows an effective reduction of AFM 
order also in presence of strong hole-phonon interactions. 



In conclusion, we demonstrated the important role 
played by the long range interactions in the problem of 
holc-phonon coupling for the physics of high temperature 
superconductors. Long-range electron hoppings lead to a 
strong anisotropy of the EPI which is possibly observed 
in the angle resolved spectroscopy of cuprates. The most 
important results are that the non-local EPI makes po- 
larons lighter and occurs to be not effective in the sup- 
pression of the motion of the hole over the NN making, 
thus, the models with strong nonlocal EPI justified for 
description of weakly doped high temperature supercon- 
ductors. 
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